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The observed prevalence of oscillatory signals in the spectroscopy of biological light-harvesting 
complexes at ambient temperatures has led to a search for mechanisms supporting coherent transport 
through larger molecules in noisy environments. We demonstrate a generic mechanism supporting 
long-lasting electronic coherence up to 0.3 ps at a temperature of 277 K. The mechanism relies on 
two properties of the spectral density: (i) a large dissipative coupling to a continuum of higher- 
frequency vibrations required for efficient transport and (ii) a small slope of the spectral density at 
zero frequency. 
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Long-lasting oscillatory signals in optically excited 
light-harvesting complexes (LHCs) point towards the 
prevalence of a coherent electronic dynamics in molec- 
ular networks at physiological temperatures [2|-|5[. The 
experiments probe the electronic dynamics using two- 
dimensional (2d) echo-spectroscopy for a series of delay- 
times between the excitation pulse and the probe pulse. 
The 2d spectroscopy makes studies of the dynamics of 
dissipative systems possible and has found further ap- 
plications in mcsoscopic systems such as molecular nan- 
otubes [6| and semiconductor devices [7|. Due to its 
known crystallographic structure and relative simplic- 
ity, the Fenna-Matthews-Olson (FMO) complex serves 
as the prototype system for studying the choreography 
of the energy transfer from the antenna to the reaction 
center of a light harvesting complex [81 . In 2d spectra 
long-lasting beatings are observed, ranging from 1.2 ps 
at T ^ 150 K to 0.3 ps at T = 277 K 0. The interplay 
of coherent dynamics, which leads to a delocalization of 
an initial excitation arriving at the FMO network from 
the antenna, and the coupling to a vibronic environment 
with slow and fast fluctuations, has lead to studies of 
environmentally assisted transport in LHCs lOl [H | . 

An important open question is whether coherence 
plays a key-role in the functioning of light-harvesting 
complexes [g . The theoretical understanding of the ex- 
periments is in its early stages and atomistic simulations 
based on molecular dynamics have not reached agree- 
ment 1^, ll3[. A calculation of 2d echo-spectra based 
on the molecular dynamics simulation [IJ] does not show 
clear coherent oscillations at T = 277 K. One key in- 
gredient for efficient transfer dynamics is the strong cou- 
plingto vibronic modes, which induces energy dissipation 
[lOl llll, ll5[ . For the FMO complex the thermalization oc- 
curs within picoseconds and was observed by Brixner et 
al. by the decay of diagonal-peak amplitudes to lower en- 
ergies [l6|. It has been proposed that the inclusion of 
the finite time scale of the reorganization process gives 
rise to long-lasting coherence in LHCs 17| . While a slug- 
gish bath relaxation leads to prolonged population beat- 



ings in the FMO network, calculations of 2d echo-spectra 
show oscillations of the exciton cross-peaks for only about 
l/6th of the experimentally recorded time at T :=; 150 K 
18| . For an even longer bath- relaxation cross-peak oscil- 
lations are absent already at T = 77 K [l9|. Hence the 
search continues for theoretical models which support a 
long-lasting coherent dynamics leading to cross-peak os- 
cillations in 2d echo-spectra and simultaneously retain 
the fast dissipation. Recently, the coupling to a super- 
position of discrete vibronic modes rather than the exis- 
tence of a purely electronic coherence has been p rop osed 
as an explanation for the detected oscillations 20, | 21 



However, no significant changes in the oscillations of the 
FMO complex have been seen in experiments which al- 
tered the vibronic modes by using native mutants, ^'^C 
substitution, and partial deuteration |22| . 

Here, we present an alternative mechanism which does 
yield long-lasting and purely electronic coherence in the 
presence of a strong dissipative coupling required for effi- 
cient transfer. The mechanism is independent of the spe- 
cific exciton-system and relies on properties of the contin- 
uum part of the vibronic spectral density as opposed to 
discrete vibronic modes. Realistic models of the exciton 
dynamics in the FMO complex have to include higher or- 
der phonon processes as well as the finite time scale of the 
reorganization process 17[. This requires to go beyond 
approximative rate equations [23[ and non-perturbative 
techniques [2J-|30J are necessary to study the dissipative 
transfer dynamics. 

One key parameter determining the duration of coher- 
ent oscillations is the spectral density, which encodes the 
the mode-dependent exciton-bath coupling. The spectral 
density has important implications for the two contribu- 
tions to the decoherence rate 7 which is the sum of the 
relaxation rate ^r and the pure-dephasing rate 7^ ([33|, 
ch. 21.4.2, m ch. 3.8.2), 



7 = 7r/2 + 7d. 



(1) 



Before we discuss decoherence-rates of the FMO com- 
plex, we analyze the impact of the spectral density on 



a two-site system with site-energies ±e/2, inter-site cou- 
plings d/2, and b"^ = e'^ + (P . The two rates in Eq. ^ 
are given by 7, « <PS{b)l{2h^), 7^ = 625(0)/(262), 
and are linked to the spectral density J{ui) via S(lo) = 
J{uj)coth{fiuj/{2kBT)). At zero frequency, S{0) is pro- 
portional to the slope of the spectral density. We eval- 
uate the decoherence rate for two different models of 
the spectral density, (i) a super-Ohmic spectral density 
Jso ^ w, and a single-peak Drude-Lorentz spectral den- 
sity Jdl(w) = IvXiojiiP' -l-o;^), depicted in the Supple- 
mentary Material (SM) [i| Fig. gH 

In the super-Ohmic case SsolO) = and the pure- 
dephasing term vanishes 7^ = and decoherence arises 
only through relaxation. The opposite regime is present 
for Jdl with 5dl(0) ~ AksTX/v. For the parame- 
ters used in models of a sluggish bath-relaxation in the 
FMO complex [17| (bath-correlation time i'~^ = 50 fs, 
reorganization energy A = 40 cni^^, e = 150 cm~^, 
d = 200 cm~^) the pure-dephasing rate dominates 
Id.DL = 0.19 cm~^K~^, over the relaxation contribution 
It.dlI^ = 0.08 cm^^K^^ to the decoherence rate. 

The spectral density for the FMO complex (circles 
in Fig. [1]) has been extracted from fluorescence-line- 
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narrowing experiments at T = 4 K [31|, |32| and shows 
a super-Ohmic behaviour. However, previous calcula- 
tions of 2d-echo spectra of seven-site FMO complex have 
been often performed with a single-peak Drude-Lorentz 
spectral density [18|, |l9[, and did not yield long-lasting 
cross-peak oscillations. The situation is different in the 
population dynamics, which has been analyzed for differ- 
ent spectral densities [26[ and yields no significant differ- 
ences in the duration of population beatings for different 
spectral density. 

As we show next, 2d echo-spectra are not directly 
derivable from population dynamics and carry along ad- 
ditional information about the exciton dynamics 35| . To 
demonstrate this, we calculate the exciton-dynamics for 
the spectral densities Jupeaks (-/speaks) shown in Fig. [l] 
which represent the experimentally determined spectral 
density by a superposition of n = 11 (n = 3) shifted 
Drude-Lorentz peaks [36[ 






(2) 



The shifts allow us to vary the slope of the spectral den- 
sity 3{y>) at w = and to mimic the super-Ohmic be- 
haviour seen in the experiment (JwendUng)- Compared 
to the single-peak Drude-Lorentz spectral-density dis- 
cussed before (A = 40 cm-\ y-^ = 50 fs, O = 0), 5(0) 
is reduced by a factor of 13 (6). Both superpositions 
■/{s.iijpeaks providc a good approximation of the spectral 
density at the transition energies corresponding to the 
exciton energy differences. 

The Hamiltonian of the FMO complex is based on the 
Frcnkcl exciton model for LHCs [341 and consists of a 
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FIG. 1: Spectral density of the FMO complex. Cir- 
cles: measured spectral density [3l|], parametrization [32 |. 
Fits J{3,ii}peaks s-re used for the GPU-HEOM calculation 
eq. ((2|. The marks on the frequency axis indicate the transi- 
tion energies corresponding to the difference of exciton fre- 
quencies. Parameters for Jspeaks: A = {10,15,13} cm~^, 
v~^ = {250,120,65} fs, Q. = {85,170,300} cm'S Jn 



{1.2, 6.4, 7.4, 15.6, 3.4, 1.8, 4, 2, 1.8, 1.9, 2} cm" 



llpcaks • 
"1 _ 



A 

{1600,550,400,370,750,800,750,600,750,750,500} fs, !^ = 

{53, 73, 117, 185, 235, 260, 285, 327, 363, 380, 480} cm-\ 



part describing the coherent exciton dynamics 

N N N 



ra=\ m n^m 



where a}„ is the creation operator of an electronic exci- 
tation of bacteriochlorophyll (BChl) m with site energy 
em- The inter-site couplings between the N BChls are 
denoted by Jmn- The BChls are coupled to vibronic 
modes modeled by a set of independent harmonic os- 
cillators Hphon = J2m,i^^M,mbi-m- The cxciton-bath 
couphng is given by Hox-phon = J^rn^'ii"'"^'^"^ ^^^^ 
Um = J2i T^^idi{bi^m + b\m). The reorganization energy 
"Hreorg = Z^m ^^m^^m with A = X]j fi^idf /2 is added to 
the site energy em- For the seven-site FMO complex 
we take the site energies and inter-site couplings from 
Ref. [32l |. For simplicity we assume identical vibronic- 
excitonic couplings at each site and neglect correlations 
of vibronic modes at different sites. The reorganization 
energy A « 40 cm~^ is of the same order as the exci- 
tonic energy-spacings, the bath-relaxation time, and the 
thermal energy ksT at T = 277 K [13, [32|. 

The strong couplings of electronic and vibronic de- 
grees require to solve the density-matrix propagation 
with a non-perturbative and non-Markovian approach. 
The computation of the 2d echo-signal for the FMO com- 
plex necessitates about 10'^ more time-steps than the 
calculation of the population dynamics and has only re- 
cently been achieved for the non-Markovian case [ISJ, ll9| . 
We solve the exciton dynamics within the hierarchical 



(a) 



CP (1.5) 



CP (1.5)-SE 



RP signal, FWHM= cm-i 

1 2 3 45 6 7 



12700 



12100 




(b) 



12100 12300 12500 12700 

Wi (CITI^^) 

SE-RP signal, FWHM= 94 cm-^ 
1 2 3 45 6 7 



12500 



3 12300 



12100 - 




12100 12300 12500 12700 



FIG. 2: (color online) Calculated 2d echo-spectra of the FMO 
complex at T = 150 K for delay time Tdoiay = 400 fs, spectral- 
density Jspcaks- (a) Real part of the rephasing signal (sum of 
ground-state bleach, stimulated emission, and excited state- 
absorption) without static disorder (linear scale). The grid of 
exciton eigenenergies is labeled 1 — 7 and cross-peak CP(1,5) 
marked by the arrow, (b) Real part of the rephasing signal 
(stimulated emission) for static disorder FWHM= 94 cm~^. 
The ArcSinh scale facilitates comparison with [0], Fig. 1. 



equation of motion (HEOM) formalism 2J, |25|. We gen- 
eralize this formalism from its original single-peak Drude- 
Lorentz term to a superposition of shifted peaks, eq. ^. 
The implementation on high- per formance graphics pro- 
cessing units (GPU-HEOM) [15[ allows us to advance 
500 coupled hierarchical equations in lock-step and to 
obtain the computationally demanding 2d echo-spectra. 
We verified the convergence of GPU-HEOM by increas- 
ing the hierarchy-depth and inclusion of low-temperature 
correction terms [171, |l8[ . 

A typical 2d FMO-spectrum at T = 150 K and delay 
time Tdeiay = 400 fs for the spectral density Jspoak is 
shown in Fig. [2ja) . The 2d echo-spectra are calculated 
as detailed in Ref. [18[ by evaluating the third-order re- 
sponse function, including excited-state absorption [37| 
for the non-disorder averaged results. We mimic rota- 
tional averaging by probing along eight molecular orien- 
tations with respect to the laser-field polarization. 
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FIG. 3: Cross-peak CP(1,5) oscillations of the FMO complex, 
real part of the rephasing signal. The sum of ground-state 
bleach, stimulated emission, and excited state- absorption is 
shown in (a) T = 150 K and (c) T = 277 K. (b) Stimulated 
emission pathway at T = 150 K for different disorder distri- 
butions FWHM= {0,47,94} cm"^ (d) Stimulated emission 
pathway at T = 277 K (no disorder). 



The inclusion of static disorder leads to a broadening 
of the peaks as shown in Fig. [IJb) for the stimulated- 
emission rephasing pathway and a Gaussian disorder- 
ensemble of 200 realizations with FWHM= 94 cm~^ 
added to the diagonal site-energies. The disordered peaks 
in Fig. 2(b) (which does not include the excited-state 
absorption) are somewhat narrower than in the exper- 
imental data 0, Fig. IC. In Fig. ^ we reference the 
peaks on the diagonal (DP) and the cross-peaks (CP) by 
their location on the exciton eigenenergy-grid. Disorder 
leads to a broad central region where the peaks coalesce 
(compare Fig.[2ja) and (b)), and amplitude is present in 
regions A, B, and C, in reasonable agreement the exper- 
imental data in [9[, Fig. IC. The calculated 2d-spectra 
differ significantly from the single-peak Drude-Lorcntz 
results shown in the SM [1[ , Fig. S5 and in [l8| . Besides 
the broader line-shapes in the single-peak Drude-Lorentz 
case, the dynamics of the diagonal-peak amplitudes is al- 
tered and an upward amplitude transfer from DP(2,2) to 
DP(3,3) is visible for Jspeaks, which is in line with the 
experimental data at T = 77 K 16[, but does not exist 
in the single-peak Drude-Lorentz case. 

To determine the coherence-time of cross-peak oscil- 
lations, we analyze the changes in the amplitude of the 
cross-peak CP(1,5) marked with the arrow in Fig. [^Ja). 
The oscillatory signal of CP(1,5) for Jspeaks shows mod- 
ulations for the rephasing pathways (Fig. [3l^a)). The 
coherent nature of the oscillations is best analyzed in 
the stimulated-emission (SE) pathway in Fig. El^b) which 
shows oscillations at T = 150 K (277 K) for delay-times 
Tdoiay up to 0.6 ps (0.3 ps). Thc excited-state absorption 
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FIG. 4: Exciton dynamics for the FMO with the spectral 
densities Jspeaks and Jupcaks- Time evolution of the coher- 
ence pEi,Ez{t) = {E\\p{t)\E5) . The duration of coherent os- 
cillations of cross-peak CP(1,5) shown in Fig. [3] is in good 
agreement with the coherence life times pEi,E^{t)- 



included in Fig. [3l[a) neglects correlation between exci- 
tons and almost compensates the stimulated-emission, 
resulting in a larger non-oscillatory background. The 
addition of static disorder diminishes the oscillations to 
some extent, but the oscillatory signal prevails as shown 
for the SE-pathways in Fig. W^) at T = 150 K. In the 
single-peak Drude-Lorentz case discussed in [ISl] and the 
SM [l[, Fig. S4, the oscillations decay two times faster 
and no long-lasting oscillations remain in the stimulated- 
emission signal at T = 277 K. 

The information of the exciton lifetime is already in- 
corporated in the simpler to calculate single exciton dy- 
namics. To this end, we prepare the system in an ini- 
tial density matrix p(t = 0) = \Ei){Ej\ given in terms 
of two eigenstates {\Ei),\Ej)^ of %cx- In the 2d echo- 
spectra setup such a delocalized exciton state is reached 
after the first two laser pules illuminate the sample. 
The life time of the oscillatory components of cross-peak 
CP(1,5) of the 2d echo-spectra in Fig. |3]is in good agree- 
ment with the life time of the corresponding coherences 
pE^^E,[t) = {Ei\p{t)\E5) shown in Fig. HI 

To isolate the impact of the vibronic mode-distribution 
we compare the coherences for the spectral densities 
>^3peaks and Jupoaks in Fig- SI The coherences are al- 
most unchanged by switching from Jspeaks to Jiipeaks- 
Specific modes in the peaked spectral density as the one 
near 180 cm~^ induce beatings in the coherences with 
small amplitude < 0.01. In contrast to (5-shaped vi- 
bronic modes, the peaks in Jupcaks which mimic the 
experimental data, have only a minor influence to the 
coherences. However, this does not imply that higher- 
frequency modes are unimportant. The relaxation to 
thermal equilibrium, which is prominently visible in 2d- 
echo spectra, relies on a large enough coupling of the 
exciton-system to the vibronic continuum. The spec- 
tral densities studied have a reorganization energy A > 
35 cm~^ which lead to a thermalization after 5 ps at 
T = 77 K as observed for the FMO complex [l^. 



For efficient energy-transfer it is important to verify 
that the decoherence rate 7 still supports a thermaliza- 
tion of the exciton system within a few picoseconds. This 
is indeed the case for Jspeaks and Jupeaks and shows that 
despite the reduced slope of the spectral density, the sys- 
tem is still strongly coupled to the vibronic environment. 
Compared to the single-peak Drude-Lorentz form, the 
spectral densities ^{3,ii}pcaks support even a faster trans- 
fer from the initial BChl 6 to the target BChl 3. In agree- 
ment with the experimental observations, the 2d spectra 
calculated here for the spectral density Jspcaks do support 
long-lasting oscillations of the cross-peaks. 

In conclusion, we have shown that two seemingly 
contradictory requirements, namely (i) strong coupling 
needed for a fast thermalization and (ii) the observation 
of long-lasting coherent oscillations are two sides of the 
same coin. The careful modeling of the continuous part 
of the spectral density towards zero frequency is essen- 
tial for the prevalence of coherence in molecular networks 
since it determines the pure-dephasing rate j^. By com- 
bining the exciton Hamiltonian with a spectral density 
with the small slope at w = seen in experiments, we 
find a much better agreement of the theoretically calcu- 
lated 2d spectra and the experimental data once static 
disorder is included. The vanishing pure-dephasing rate 
of a super-Ohmic spectral density implies that the deco- 
herence is driven by the relaxation rate 7^. This is ex- 
actly the opposite behaviour compared to models based 
on the single-peak Drude-Lorentz spectral density jl7| . 
Furthermore, we find only a very small influence of the 
additional peaks in the spectral density Jupcaks on the 
coherences. The cross-peak oscillations require a calcu- 
lation of 2d-echo spectra beyond the weak-coupling ex- 
pansion for structured spectral densities, which becomes 
feasible by the efficient GPU-HEOM algorithm. 
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Appendix 



Comparison of exciton dynamics for different spectral densities 



The exciton dynamics is strongly affected by the spectral density and the dephasing rate in Eq. (1) depends on 
the slope of the spectral density towards zero frequency. Fig. [5] shows a close up of spectral densities discussed in 
the main text, and in additions includes for convenient reference the single-peak Drudc-Lorcntz density (dotted line, 
>/dl), used in some previous studies |17H19l |. 
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FIG. 5: Left panel: low frequency onset of the spectral density of the FMO complex. Right panel: global view. Circles: 
measured spectral density [3l||, parametrization [S^l. Solid line Jspcaks, dashed line Jupcaks, dotted line Jdl, parameters given 
in the main text. 
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FIG. 6: Detail of the long-time evolution of the coherence pEi,Er^{t) ~ {Ei\p{t)\E5), which shows small-amplitude oscillations 
for Jiipcaks related to the peaks of the spectral density Jupcaks located at higher frequencies. For Jdl (dotted line), almost no 
coherences prevail. 



Population dynamics and coherences 



The additional peaks present in Jupoaks at higher frequencies, which are not contained in Jspcaks, affect the coher- 
ences of the off-diagonal elements pei,e^ only weakly, Fig. [6] The influence of the additional peaks in the spectral 
density Jupoaks becomes visible as small-amplitude oscillations for longer delay times (Fig. [S]). This shows how the 
initially electronic coherence gets transformed into oscillations related to the details of the sharp peaks in the spectral 
density Jupoaks- This finding is in-line with the theoretical calculations for spectral-densities with added J-peaked 
vibronic modes [20|. The same observation applies to the relaxation process, which is shown in Fig.[7]for a population 
dynamics starting either at site 1 or at site 6. For comparison, the results for the Drude-Lorentz spectral density 
(dotted curve) is also included, which deviates most significantly for pge- 
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FIG. 7: Population dynamics for an initial population at sites 1 and 6 for the 3 spectral densities shown in Fig. [5] and for 
temperatures T = 150 K and T = 277 K. The relaxation proceeds on very-similar time-scales for the two spectral densities 
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FIG. 8: Exciton dynamics for the FMO complex for two different spectral densities Jspcaks and Jdl. (a) Real part of cross- 
peak CP(1,5) of the signal (rephasing and non-rephasing) at T = 150 K and T = 277 K. (b) Coherent oscillations of the 
stimulated emission (SE) part of cross-peak CP(1,5). 

Differences in 2d spectra 

The cross-peak oscillations depend strongly on the initial slope of the spectral densities towards zero frequency. 
This is reflected in Fig. [SJ which compares the cross-peak oscillations for the two different spectral densities Jspcaks 
and Jdl, for which the dephasing- and relaxation contributions to the decoherence rate are reversed in magnitude. 
Fig. |Hl[b) shows that for Jdl (dotted lines) the cross-peak oscillations are largely reduced at T = 150 K and almost 
completely absent at T = 277 K. while the super-Ohmic onset of Jspcaks still leads to visible oscillations at T = 277 K 
(solid lines). Moreover, also the amplitudes of the diagonal-peaks show pronounced differences for the two spectral 
densities. The spectral density Jspcaks transfers the highest diagonal-peak amplitude upwards with increasing delay 
time from DP(2,2) to DP(3,3) as illustrated in the left panels of Fig. |8l This shift is also seen in the experimental 
data at T = 77 K [ifll, Figure 1, panels a,b. 

In the Drude-Lorentz case (Fig. |51 right panels) the highest amplitude is located at diagonal-peak DP(2,2) for all 
delay-times, which corresponds to previous models (see Fig. 2 in j38l|). but differs from the experimental data and the 
result for Jspeaks- 
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FIG. 9: (color online) Calculated 2d echo-spectrum of the FMO complex at T = 277 K for the delay times Tdoiay = {40, 100} fs. 
Shown is the real part of the rephasing and non-rephasing signal. The square marks the diagonal peak (DP) with the highest 
amplitude. Left panels: spectral density Jupeaks, showing amplitude shifts from DP(2,2) at Tdday = 40 fs to DP(3,3) for 
Jdeiay = 100 fs, right pauels: spectral density Jdl, highest amplitude remains at DP(2,2). 



